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ABSTRACT 

We report results of a series of non radiative N-body/SPH simulations in a KCDM cosmology, designed to 
study the growth of angular momentum in galaxy systems. A sample of 41 halos of differing mass and envi- 
ronment were selected from a cosmological N-body simulation of size 32.5/!~'Mpc, and re-simulated at higher 
resolution with the tree-SPH code GADGET. 

We find that the spin of the baryonic component correlates well with the spin of the dark matter, but there is 
a misalignment of typically 20° between these two components. The spin of the baryonic component is also 
on average larger than that of the dark matter component and we find this effect to be more pronounced at 
lower redshifts. A significant fraction / of gas has negative angular momentum and this fraction is found to 
increase with redshift. This trend can be explained as a result of increasing thermalization of the virializing gas 
with decreasing redshift. We describe a toy model in which the tangential velocities of particles are smeared 
by Gaussian random motions. This model is successful in explaining some of the global angular momentum 
properties, in particular the anti-correlation of / with the spin parameter A , and the shape of the angular 
momentum distributions. 

We investigate in detail the angular momentum distributions (AMDs) of the gas and the dark matter components 
of the halo. We compare and contrast various techniques to determine the AMDs. We show that broadening 
of velocity dispersions is unsuitable for making comparisons between gas and dark matter AMDs because 
the shape of the broadened AMDs is predominantly determined by the dispersion and is insensitive to the 
underlying non broadened AMD. In order to bring both gas and dark matter to the same footing, we smooth 
the angular momentum of the particles over a fixed number of neighbors. The AMDs obtained by this method 
have a smooth and extended truncation as compared to earlier methods. We find that an analytical function in 
which the differential distribution of specific angular momentum j is given by P{j) = .^.L A j)"'^^'^^''''^ where 

id = jtotjoL, can be used to describe a wide variety of profiles, with just one parameter a. The distribution of 
the shape parameter a for both gas and dark matter follows roughly a log-normal distribution. The mean and 
standard deviation of log(a) for gas is -0.04 and 0.11 respectively. About 90-95% of halos have a < 1 .3, while 
exponential disks in NEW halos would require 1.3 < a < 1.6. This implies that a typical halo in simulations 
has an excess of low angular momentum material as compared to that of observed exponential disks, a result 
which is consistent with the findings of earlier works, a for gas is correlated with that for dark matter (DM) but 
they have a significant scatter < acas/ctOM >= 1 -09 ± 0.2. acas is also biased towards slightly higher values 
compared to odm- The angular momentum in halos is also found to have a significant spatial asymmetry with 
the asymmetry being more pronounced for dark matter 
Subject headings: cosmology: dark matter — galaxies: formation — galaxies: structure 



1. INTRODUCTION 

Disk galaxies are rotationally supported systems and their 
structural properties are intimately linked to their angular mo- 
mentum distribution. The standard picture of formation of 
disk galaxies is that the density perturbations grow due to 
gravitational instability and end up forming virialized sys- 
tems of dark matter and gas. The gas cools and collapses 
towards the center (White & Rees 1978) . The gas has angular 
momentum which it acquires due to tidal interactions. This 
can be quantified in terms of a dimensionless spin param- 
eter A = J\E\^ ^^/(GM^^^) (Peeble s 1969 ) which has a value 
of about 0.05 JEfstathiou & Jonesll979tlBarnes & Efstathioul 
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ll987HSteinmetz & Bartelmann I I 19951) . The angular momen- 
tum of the gas is conserved during the collapse resulting in the 
formation of a centrifugally supported disk, whose size is con- 
sisten t with that of observed disk galaxies (iFall & Efstathioul 
fT98(il) . 

Based on the initial density profiles (p = p{r)) and angu- 
lar momentum distributions (m = m(y)) or (M<j = M^jij)) 
the final surface density of disks can be determined. 
Dalca nton. Soergel. & Summers (1997) derived the surface 
density of the disks assuming the halos to be uniform spheres 
in solid body rotation. From numerical simulations it is now 
known that density profiles of dark matter in halos follow 
a universal p rofile as described by Navarro. Frenk. & Whit^ 
(119961 11997 . NEW). Using these realistic profiles and as- 
suming the final surface d ensity of disks to be exponential 
iMo. Mao. & Whita (119981) investigated various properties of 



disks like rotation curves, disk scale lengths and so on. They 
also addressed the issue of stability of disks. However, an- 
gular momentum distributions are required to more realisti- 
cally m odel the surface density of disk galaxies. Bullock et al. 
(|200r, henceforth B2001) have reported that in CDM simu- 
lations, the DM halos obey a universal angular momentum 
distribution of the form 



M{< j) nj 



M.. 



Jo + J 



(1) 



where j is the specific angular momentum and M(< j) is mass 
with specific angular momentum less than /'. The shape pa- 
rameter n has a log-normal distribution and the 90% range is 
given by 1 .05 < /i < 2.0. Assuming that the angular momen- 
tum profile of gas is identical to that of dark matter, B2001 
calculated the surface density profiles of the resulting disks, 
and found that (for the range of ^ given above) the result- 
ing disks are too centrally concentrated compared to expo- 
nential profiles. In addition, detailed hydrodynamical simula- 
tions of gas collapse in hierarchical structure formation sce- 
nario exhibit the so-called "angular momentum catastrophe": 
the gas component looses its angular momentum due to dy- 
namical friction and ends up forming disks that are far too 
concentrated (Navarro & Benz* 1991*: 'Navarro & Whit Jl994t 
JNavarro & Steinmetz 1997; Steinmetz & Navarro 199^. 

Recently van den Bosch et al. (2002) (henceforth vB2002) 
has tested the assumption, that the AMDs of gas and DM 
are similar, by directly measuring the AMDs of gas in proto 
galaxies from hydrodynamical simulations at z = 3. They find 
the spin parameter distribution of gas and dark matter to be 
identical in spite of the angular momentum vectors of gas and 
dark matter being misaligned by ^^ 35°. In order to compare 
the AMDs of gas and DM within halos, they broaden the ve- 
locities of gas to account for the microscopic random motions 
of gas atoms. The broadened profiles are found to be very 
similar to that of DM, though, as we are arguing later in this 
paper, this result is dominated by the dispersion hiding away 
the effect of the actual profiles. 

vB2002 also demonstrated that a considerable amount of 
gas (between 5 and 50 percent) have negative angular mo- 
mentum. Assuming that the gas with negative angular mo- 
mentum combines with that of positive angular momentum to 
form a non rotating bulge, and the remaining positive angular 
momentum material ends up forming a disk, they found the 
surface density profiles of the resulting disks to closely fol- 
low an exponential profile. However the galaxies end up with 
a large B/D ratio, and with a minimum B/D of 0. 1 the ques- 
tion of how to form bulge-less dwarf and LSB galaxies is still 
unanswered. The problem of an excess of low AM material 
gets transferred into a problem of excessive bulge formation. 

To investigate these issues in more detail and to see if these 
properties have any evolution with redshift we perform a se- 
ries of N-Body/SPH simulations of selected halos till z=0. 
Special care is taken to have high number of particles in the 
final virialized halos so as to measure the angular momentum 
accurately. The details of the simulation and methods of anal- 
ysis are described in Sec-2 . Results related to global angular 
momentum parameters and their evolution with redshift are 
presented in Sec-3. In Sec-4 we present a toy model to ex- 
plain some of these findings. In Sec-5 and 6 we analyze the 
angular momentum distributions. 

2. METHOD 



2.1. Simulation 

The cosmology adopted in the simulation is the so called 
concordance model of ACDM cosmology, in agreement 
with recent WMAP and SDSS results ( Soergel et al. 2003t 
iMelchiorri. B ode. Ba hcall. & Silk! ^2003^ ■ The parameters 
adopted are fl\ = 0.7 ,n„, = 0.3 and 17^ = 0.02235/?"^ and a 
Hubble constant of 65kms"' Mpc~' . The power spectrum 
parameter determining the amplitude of mass fluctuations in 
a sphere of radius 8/i~' Mpc was set to 0.9 and shape parame- 
ter r was set to 0.2. An AP^M code was used to evolve 128^ 
dark matter particles in a 32.5/i~' Mpc cube from z = 24 to z = 
using 2000 equal steps in expansion factor At z = 0, 41 halos 
were selected with circular velocities ranging from 64kms"' 
to 310 kms"' and masses ranging from 1.3 x 10" M0 to 
1.5 X IO^-'Mq . These halos were then re-simulated at higher 
resolution with more dark matter particles and also with an 
equal number of gas particles. The simulations were per- 
formed from z = 50 to z = 0. T hese re-simulat ions were done 
using the code GADGET ( Soringel. Yoshida. & White-20011) . 
The gas particles were given an initial temperature of \00K 
at z = 50 and an artificial temperature floor of lOOK was kept 
during the simulation. The number of particles (of each kind 
) within the virial radius ranges from 8000-80000. A gravi- 
tational softening of 2kpc /i"' (physical) was used. The inte- 
gration was performed in comoving co-ordinates. 

2.2. Halo Identification 

We adopt the method of vB2002 to identify the center of 
mass of a virialized region. We start with the densest gas 
particle as a guess for the center of mass and iteratively in- 
crease the radius till the average mass density enclosed by 
a spherical region is A,, times the mean matter density at 
that redshift. A,, is approximated by ("Br van & Normanll998h 
A,, ~ {I^t:^ + i2x + -1,9x^)1 {I +x),whtx&x=n„,(z)-\. We 
re-center the particles within the virial radius in velocities and 
position and then recalculate the virial radius based on this 
new center of mass. We repeat this process until the distance 
between the center of mass before and after the calculation of 
virial radius is less than 0. 1 percent of the virial radius. 

2.3. Angular momentum distributions 

As discussed by vB2002, there are two kinds of velocities 
for particles in the simulations, the actual microscopic veloc- 
ity V of individual particles and the mean streaming velocity 
u at any location x . The actual microscopic velocity v is 
given by equation v = u + w, where w is the particles random 
motion. For collisionless dark matter particles, which interact 
only through gravity, the velocity given by simulations is v, 
whereas coUisional gas (SPH) particles, the velocity given by 
simulations is u, the information about the random motion is 
incorporated into the internal energy per unit mass U . If a is 
one dimensional velocity dispersion of the particle,its temper- 
ature T is given by 

3 , l>kT 

n'^ -n — 

where /i is the mean molecular weight of gas. In order to com- 
pare the kinematical properties of the gaseous or dark matter 
component we either need to broaden the velocities of gas 
particles by using Eq. (|3 (we label these by superscript t; de- 
noting the actual motion ) or smoothen the velocities of dark 
matter particles by an appropriate smoothening length (we la- 
bel these by superscript s; denoting the streaming motion of 
flie fluid). 



U=-a' 



(2) 



The total angular momentum of gas or dark matter is given 



by 



Kiis.Ndm 



I gas, DM '■ 



y^ OT,ri X Vi 



(3) 



/=! 



where ri and y-, are the radius and velocity vectors respec- 
tively of a particle /, in a co-ordinate system in which both the 
position and the velocity of the center of mass of the entire 
halo (DM + gas) is zero. For the spin parameter, we use the 
modified definition of B2001 . 



A 



'gas, DM — 



\jgas,DM\ 
V2R,,,.V,ir 



(4) 



where jg„s and jdm are the mean specific angular mo- 
mentum of gas and dark matter respectively, and y,,„- = 
^G{Mgas + MDM)/Rvir IS the circular velocity at virial radius 
Rvir- The misalignment between the angular momentum 
vectors of gas and dark matter Jgas and Jdm is given by 



■■{■ 



I gas 



Jdm 



Jgas 1 1 Jd: 



(5) 
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We make a co-ordinate transformation such that the z-axis 
is aligned with the total angular momentum vector For gas 
particles the z-axis is aligned along Jgas while for dark mat- 
ter z-axis is aligned along Jdm- The z component of specific 
angular momentum jy is then measured (henceforth we drop 
the subscript and denote it by f). The fraction of mass with 
y < is labeled as fgas and foM for gas and dark matter re- 
spectively. The differential angular momentum distribution 
(AMD) is the fraction of mass P{j) with specific angular mo- 
mentum between j to j + dj i.e. j_ P(j)dj = 1. We define a 
parameter / which is related to j by 



I: 



J 



V2VyirR.i 



(6) 



(similar to definition in vB2002 except for a factor of V^). 
The above definition for I implies that 



lPil)dl = A 



(7) 



The total specific angular momentum of a halo is denoted 
by jioi and this is related to A by j,o, = \V2RrirVyir, see 
Eq. 0. The shape and ex tent of profiles depend on it„t . With 
s = i I i ini ,as proposed in Ivan den Bosch. Burkert. & SwatersI 
J2001L henceforth BBSOl), we obtain 



sP(s)ds = 1 



(8) 



In most of our analysis we neglect the negative tail and de- 
fine the distributions for the positive tail only and normalize 
with respect to it. The range of s or /' in that case is from to 
oo and for halos with negative tails the quantities like j,oi also 
need to be recalculated for the positive tail only. The cumula- 
tive angular momentum distribution P{< j) is the fraction of 
mass with AM less than equal to j. It is also defined only for 
the positive tail, and is normalized with respect to it. 

3. RESULTS 

3.1. Global angular momentum parameters at z=3 

As a first test we compare the results of our simulations 
against the analysis of vB2002, who compared various global 
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Fig. 1 . — A comparison of various global angular momentum parameters 
at z = 3 for both gas and dark matter. The parameters compared are the spin 
parameter A and the fraction with negative angular momentum /. A compar- 
ison with Fig-2 of vB2002 shows that at z = 3 the results of our simulations 
are in good agreement with the analysis of vB2002. 



angular momentum parameters using a sample of 378 halos 
at a redshift of z=3 ( see Fig-2 in vB2002). Our results for a 
sample of 41 halos at the same redshift, are shown in Fig.^ 
We find them to be in good agreement. 

The first plot on upper left compares Xgas against Xdm- A of 
gas and DM are well correlated with a Spearman rank coef- 
ficient of r, = 0.88, but show significant scatter The mean of 
^gas/^DM is 1-27 with a standard deviation of 0.40. This sug- 
gests that statistically the angular momentum acquired by gas 
and dark matter is similar but on a one-by-one comparison the 
spin parameters can be quite different. The fraction of mat- 
ter with negative angular momentum for gas and dark matter 
denoted by f and f^/^, are also well correlated with Spear- 
man rank coefficient of 0.81 and < /'dm/ fgas >= 0-96 ± 0.07 
(As described in Sec-2.3 the superscript t denotes inclusion 
of microscopic random motions while superscript s stands for 
streaming motions only). In the middle panels the strong anti- 
correlation of /^^j and f^j^ with their respective spin parame- 
ters can be seen. In the lower left panel /* ^ is also found to be 
anti-correlated with Xgas- If we assume that the random mo- 
tions are responsible for the negative angular momentum and 
that the ordered motion contributes to the A, then this result 
is easy to understand. The more the ordered motion the less 
will be the effect of random motions. The plot of /* ^ vs Xgas 
is found to have more scatter than that of f' vs Xgas which 
has broadened velocities. We will follow up on this though in 
more detail in section 4. 

In the lower left plot /^^^ is shown against Z^^,. /^^^ is less 
than f' as expected ( because f^as has more random motion) 
and the difference /* ^ - f is found to increase at lower / 
values or at lower A, respectively. A comparison of some of 
our results with vB2002 is given in Table [J 



3.2. Redshift dependence of angular momentum parameters 



Table 1 . Global angular momentum parameters 
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< \as - >^DM > 


-0.001 ±0.020 


0.006 ±0.009 


< fgas/lDM > 


0.97±0.10 


0.96 ±0.07 


<e> 


36.2 


18.9 



Note. — The properties shown above ai'e for halos at z = 3, alongside are 
results from vB2002 also at same redshift. The results reported in this paper 
are in good agreement with vB2002 
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Fig. 2. — Distribution of spin parameter A for dark matter (left) and gas 
(right) for a sample of 41 halos for redshifts ranging from 4 to 0. The solid 
line isthe best lit log-normal distribution as given by Eq. |5). The fit param- 
eters A and <tx are also labeled on each of the plots. 



Distribution of Xdm and Xgas 

It has been suggested by numerous studies that the distribu- 
tion of A can be described by a log-normal distribution 



p{X)dX = 



1 



Inax 



exp 



ln2(A/A) 



2al 



dX 

A 



^- (9) 



Ivan den Bosch et all (l2002h have found that the gas and dark 
matter have very similardistribution of spin parameters. They 



found A„ 



: 0.039 and Ar 



: 0.040. 



In Fig.|2]we have plotted the distribution of A^^j and Xom 
for redshifts from to 4 . It is observed that Xom fluctuates 
between 0.029 to 0.034 while Xgas seems to increase from 
0.034 to 0.041 with decrease in redshift. The fact that A^,,, is 
higher than Xqm can be seen more clearly in Fig.|3]where A^o, 
is plotted against Xdm- < Xgas/ Xdm > is greater than 1 and 
increases as z decreases. 
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Fig. 3. — Ag„j vs Xom for various redshifts. < Xgas/XoM > and < Ajaj- 
AflM > along with their standard deviations are shown on each plot. 
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Distribution of Misalignment Angle 9 

In Fig. 13 we show the misalignment angle between the total 
angular momentum vectors of dark matter and gas for various 
redshifts. The mean and median values are also shown in each 



Fig. 4. — Distribution of misaUgnment angle 9 for various redshifts. The 
dashed line is the expected distiibution if the angle is completely random or 
uncorrected. 
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Fig. 5. — Counter-rotating fraction / as a function of spin parameter A for various redsliifts. Superscript t stands for analysis done witli velocit ies w liicli 
incorporate random motions and .s stands for analysis done with streaming velocities. The solid lines are the fits given by the functional form in Eq. <10l . The 
best value of parameter Aq is indicated on each plot. 



panel. The distribution does not show any significant change 
with redshift. The mean value < > is around 20 degrees. 

Correlation of f with A 

As mentioned in Fig.^ the fraction of matter with negative 
angular momentum / is anti-correlated with A. We found that 
the / vs A distribution can be well fit by a function of the form 



A 






(10) 



/Ittct 



■dt is a Gaussian integral. The results 



where Ig(x) = J^^ 

of the applied fit are shown on each plot in Fig.|5] The anti- 
correlation can be described by a single parameter Aq. For 
the actual motion of both gas and DM, it is found that Ao is 
around 0.17 and does not seem to change much except for a 
decrease in its value at z=0. For streaming motion of gas the 
parameter Ao decreases consistently with decrease in redshift, 
thus hinting to a systematical change in the distribution of 
fgas with redshift. Furthermore, the scatter in / vs A plots 
corresponding to streaming motion of gas ( column 3 Fig.|5} 
is quite large specially at lower redshifts as compared to the 
plots for broadened motions (column 1 and 2 Fig.|5}. 

3.3. Variation of < fgas > with thermalization of gas 

In Fig.|6lthe distribution of/ is plotted for various redshifts. 
It can be seen that the distribution for dark matter does not 
seem to change with redshift, while for gas it shifts towards 
smaller fractions of negative angular momentum material at 
lower redshifts. For DM < / >~ 0.40 while for gas < / > 
decreases from 0.33 to 0.15 monotonically with decrease in 
redshift (Fig. 0. The kinetic energy of gas particles in SPH 
simulation, is a combination of translational energy and in- 
ternal (thermal) energy. The ratio Fjr = KEjr / {KETr + KEjh) 
gives an estimate of thermalization. The lower the Fjr the 
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Fig. 6. — Distribution of counter-rotating fraction / for gas and dark matter 
for a sample of 41 halos. The distribution is almost independent of redshift 
for dark matter while for gas / is lower at lower redshifts. The solid line is 
the distribution as predicted by Eq. 1251 . 



greater the thermalization. Fig.0shows that < Fjr > just like 
< fgas >, also decreases monotonically with decrease of red- 
shift. The correlation of < Frr > with < fgas > is easy to 
understand. 

The dispersion in velocity of particles gives rise to particles 
with, negative angular momentum. For very small velocity 
dispersion of particles the fraction with negative angular mo- 
mentum f ^ 0. If velocity dispersion is very large or AM is 
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Fig. 7. — Counter-rotating fractions < fgas > , < fdm > and tlie fraction of 
energy for gas in the form of translational motion < Ft,- >, as a function of 
redshift. The average is taken over a sample of 41 halos. At lower redshifts 
the gas gets more and more thermalized, its kinetic energy getting converted 
into thermal, and this decreases < / >. 



very small, then / ^ 0.5. The SPH particles only have macro- 
scopic flow velocities, the thermal energy is incorporated into 
U. For DM particles all the energy is in the form of veloci- 
ties of particles. So fgas is always less than fdm- At higher 
redshift there are more mergers and the gas is more turbulent. 
As the redshift decreases, the gas undergoes relaxation and 
kinetic energy of gas gets converted into thermal energy. The 
velocity dispersion of gas decreases resulting in a decrement 
of /g„5. For dark matter Ft,- = 1 and < / > does not show any 
considerable evolution. 

3.4. Effect of numerical resolution 

To check whether the numerical resolution has an effect on 
fgas we plot fgas vcrsus A^, the number of gas particles in the 
halo (filled squares in Fig.|8} and find no correlation. To check 
the effect of numerical resolution in more detail 30 halos were 
re-simulated at a lower resolution with 1 /8 times the original 
number of particles (open triangles in Fig.|8}. The values of/ 
do not show any systematic trend with change of resolution. 
Halos simulated in lower resolution have distributions of A ,/ 
and misalignment angle 6 (Fig.|9} identical to the distributions 
found in high resolution simulations. So the results shown 
above are robust to the effect of numerical resolution. 

3.5. Evolution of gas particles having negative angular 
momentum 

We investigate whether the gas that ends up with negative 
angular momentum also had negative angular momentum in 
the past. We identify a region in past that ends up in the virial- 
ized halo at z=0 by simply tracking halo particles at z=0 back 
in time. We re-center them, and then calculate their angular 
momentum. The fraction of particles with negative AM at any 
given stage during the evolution of this Lagrangian volume is 
a function of redshift so we denote it by fgiohaiia)- We iden- 
tify a subset of particles that were counter rotating at z = 0, 
and denote them by subscript z = 0. The fraction of particles 
out of this subset that are counter rotating at any given in- 
stant is denoted by f~o{a). Similarly /_i5(fl) is the fraction 
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times the original resolution also do not show any systematic shift. 
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Fig. 9. — Various properties of halos and their distributions at z = from 
low resolution simulations (1/8 times the normal resolution). A comparison 
with Fig.|2] Fig.|3] Fig.|4], Fig.l6lwhich are the con'esponding plots with high 
resolution simulations shows no difference. 



of counter rotating particles out of a subset of particles that 
were counter rotating at z=15. If the particles with negative 
angular momentum at z = were also counter rotating at an 
earlier epoch, then they would have /,=o(fl) ^ 1 independent 
of redshift. Similarly if counter rotating particles atz= 15 also 
had negative angular momentum later during their evolution 
then fz=i5{a) ^ 1. On the other hand if /-.=i5 or /,=o ^ f global 
then the particles are just a random subset drawn from the 
original halo and have no relation with their past or future, 
respectively. 




0.4 0.6 0.8 1.0 

Scale Factor a 
Fig. 10. — The evolution of fraction / of particles with negative angular 
momentum with scale factor a for various sub-components of gas. Solid line 

- /global vs "> Dotted line(with stars) - /_o ^s " and Dashed line(with squares) 

- /j=i5 vs a. /global is the fraction of gas with negative angular momentum for 
the whole halo, /--o is the coiTesponding fraction out of a subset of particles 
that had negative angular momentum at z = 0. /;=i5 is the fraction out of a 
subset of particles that had negative angular momentum at z = 15. 



In Fig.[ro|we have plotted fgiohai, fz=o and/— 15 as a function 
of scale factor a of the universe for a randomly selected halo. 
/j=o and /j.=i5 are both close to 1 at their respective ends but 
by z = 1 both fractions have dropped to f global and continue to 
remain so in respective directions. This suggests that particles 
that are counter rotating at a given instant will not necessarily 
counter rotate at a later epoch in future but will instead get 
mixed up randomly with the remaining portion of the halo. 

4. A TOY MODEL: GAUSSIAN SMEARING OF ORDERED 
VELOCITIES (GSOV) 

To get a better intuitive understanding of some of the results 
presented so far, in particular the anti-correlation of/ with A, 
we present here a toy model that describes the motion of the 
particles in the halo. We take a co-ordinate system with the 
z axis pointing along the direction of the total angular mo- 
mentum vector The actual velocity u of a particle consists 
of the ordered motion Vo4> (which is motion in circular orbits 
at a constant speed vv), superimposed by random motions y^. 
Each component of random motion Va is drawn out from a 
Gaussian distribution with a dispersion a. The random mo- 
tions are assumed to be isotropic, u can be written as 



U = Vo0+ VctJ+ Vo-J + Vc,„k 



(11) 



According to this model the histogram of u^ in a halo is given 
by 



P(u^)du^ 



1 



Inatf, 



-(u^-.,,f/(2al)^ 



(12) 



In Fig. E] we have plotted the histogram of u^ for various 
halos, for both dark matter and gas. For most halos the u^ dis- 
tribution of gas particles can be well described by a Gaussian. 
Some halos are biased either to right or left, and some show 
two peaks which can described by two Gaussian functions. A 
closer examination reveals that these peculiarities are associ- 
ated with substructure of halos. For dark matter a is large and 
this washes out any peculiarities that may have been present, 
consequently no significant deviation from the single Gaus- 
sian structure can be observed. Sometimes a slightly sharper 
peak compared to a Gaussian distribution can be seen indica- 
tive of a velocity dispersion that is not strictly isothermal. 




Fig. 1 1 . — A plot showing distribution of 11^ for gas and DM for four 
different halos. Solid line: the histogram for particles in simulation. Dashed 
line: A Gaussian fit to the above histogram . The shaded region corresponds 
to particles with negative angular momentum. 



In terms of this model the motion within a halo can be de- 
scribed by two parameters Vo and ( where ( = ^. C is a 
measure of random motion relative to ordered motion. We 
calculate Vg and a^ for simulated halos by fitting a Gaussian 
profile as given by Eq. J12> to the histogram of u^. Fig. ^] 
shows the distribution of C and Vo/V,, for both gas and DM at 
z = 0. logiQ can be fit by a Gaussian both for DM and gas. 
For gas C ^ 1 while for DM C ^ 5. 

If the above model of Gaussian smearing is a realistic rep- 
resentation then the fraction of matter with negative angular 
momentum is simply the probability that m^ is less than zero, 
i.e. vv^ is less than -Vo (Fig.fTst. This can be expressed by an 
integral 



/(Vo,cr0) = 



-v„ „-..V(24) 



InaA 



-dv 



= 0.5[l-£r/(^^)] 

V2cr0 



(13) 



A is related to < j^ > by 



where L(x) = f '^ , — — dt is a Gaussian integral. 

< iz > 



A: 



V2V,R, 



(14) 



If we further make the approximation that < j,/r >=< / > 
/ < r >, then using Eq. J14> we can write 

Vo = < M^ > 

= < jz/r > 



Gos solid DM-doshed 

( -0,01, 0,15) ( 0,70, 0,20) 
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Fig. 12. — Distribution of parameter C = T" ^nd -^ for a sample of 
41 halos at z=0. Distributions are roughly log-normal. Gaussian fits are 
applied and the fit parameters are indicated in upper left and right comers. 
The position of the peak in normal units is shown below the fit parameters. 



V2X/(VrRr) 



<r> 

V2XV„ 



(15) 



< r/R, > 
This can be used to express / in terms of A as shown below 

V2\ 



f=^-hi 



< r/R, > ia^/V,) 



K.)-/C(j 



(16) 



where k,- =< r/R^, > and fco- = {(T^/Vy) are defined to be two 
quantities that are constant for a halo. For an NFW halo with 
concentration parameter c 

kr = < r/Ry > 
J rp(r)d^r 
" Jp(r)d'r 

_ 7r[l+c-l/(l+c)-21n(l+c)] 
~ 4 c[l/(l+c)-l+ln(l+c)] 

~0.32 fore =10.0 



We define 



then 



Xo=<kr Xka > 



^-'-^^'Yo' 



(17) 
(18) 

(19) 



Anti-correlation off with A 

Eq. ( I19> explains the anti-correlation of / with A as dis- 
cussed in Section 3.2.3. Rightmost column in Table Hand 
Table |3llists the values of parameter Aq obtained by fitting this 
equation to the data from simulated halos ( Fig.|3. Ao as pre- 
dicted by the toy model can be calculated by using Eq. ( I18> . 
For this first kr and A;^. are calculated for each of the halos 
in the simulation then at any given redshift < kr >sim and 
< ka >sim are calculated by taking the mean over all the halos 
at that redshift. These values are listed in columns 2 and 4 
of the aforementioned tables. Ao calculated from these val- 
ues is listed in column 6. The values of Aq predicted by the 
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Fig. 13. — Dispersion of ordered velocity Vo by a Gaussian with width 



(7. P(v, a) -- 



vlTva 



-'V(2Tp 



is the Gaussian distribution function. A parti- 



cle has negative angular momentum if Va-, < —Vo, this is represented by the 
region under the Gaussian curve which is not shaded. 



Table 2. *:,- , k^ and Aq for Gas 



<kr > 



<ka > 



Ao 



Sim 



ENS 



Sim ENS pred 



Sim 



0.308 0.315 0.163 

1 0.343 0.349 0.181 

2 0.371 0.373 0.210 

3 0.380 0.390 0.217 

4 0.395 0.402 0.213 



0.036 0.036 

0.044 0.053 

0.055 0.063 

0.058 0.074 

0.059 0.074 



Table 3. kr ,ka and Aq for Dark Matter 



<kr > 



<k„ > 



Ao 



Sim ENS Sim ENS Pred Sim 






0.308 


0.315 


0.632 


0.654 


0.138 


0.145 


1 


0.350 


0.349 


0.613 


0.605 


0.152 


0.170 


2 


0.382 


0.373 


0.595 


0.578 


0.161 


0.178 


3 


0.398 


0.390 


0.590 


0.561 


0.166 


0.177 


4 


0.409 


0.402 


0.573 


0.551 


0.166 


0.189 



Note . — S im: simulations, ENS: calculated theoretically by using 
Eq. 1171 and Eq. 1211 . c used in calculations is estimated by algorithm given in 
Eke. Na van'o. & Steinmetz i2001) Ao Pred: Ao as predicted by the toy model 
Eq. 1181 . the values of < kr > and < ka- > used are the ones obtained from 
simulations, Ao Sim: An as obtained by fitting Eq. 119t to the / vs A data 
from simulations (Fig.lsl. 



toy model have a small offset but otherwise they are in agree- 
ment with those obtained from simulations. The increase of 
Ao with redshift can be understood in terms of variation of 
<kr> and < k^ > with redshift. k,- as given by Eq. ( I17> is a 
monotonically decreasing function of c. c on the other hand 
has a dependence on mass, redshift and cosmology given by 



C = C£A-5(cr8 , r , r^A , r^o , Z, M/ , Mo) 



(20) 



which can be calculated by means of an algorithm given in 
lEke. Navarr o. & Steinmetz (2001) . For calculating < c > 
at a given redshift z we put in Eq. ( I20> Mq =< M,. >=the 
mean mass of a halo at that redshift. The predicted ENS 
values of < A:r > are shown in the table, they are in good 
agreement with those obtained from simulations. For the case 
here < c >^ 1 2.0/(1 +z) (in agreement with redshift depen- 
dence given in iBullock et alJ (t2001.) ). This is the cause for 
increase of < A;^ > with redshift both for gas and DM. k^ 
for DM also depends on c because for DM the one dimen- 
sional velocity dispersion is given by cr^ = (l/3)V^fc, where 
fc ^ 2/3 + (c/21.5)"'' ( assuming that t he velocity distribu- 
tion is isotropic and homogeneous, see iMo. Mao. & Whita 

dms)). 



ka = i^ = VfJ^' 



0.65 



fore =10.0 



(21) 



The slight decrease of k^ with redshifts (Table O is again due 
to increase of < c > with redshifts. For DM kr decreases with 
c while ka- increases this makes Ao nearly a constant for all the 
halos at a given redshift and this is one of the reason for the 
small scatter in / vs A plots for DM (column 3 Fig.|5}- For gas 
a^ is not related to V,, and this results in a large scatter seen 
in its / vs A plots (column 3 Fig. |5j- But ka- can be written 
in terms of Fj,- as < k^ >~ \jFjrl'i and < Fjr > increases 
with redshift (Section 3.3 ,Fig.E}, explaining the increase of 
kcr with redshift shown in Table |2l 

It was found in B2001 that standard deviation aj of angular 
momentum for a subsample of A^ particles for DM scales like 






= Jl/N+l/{25X^N)'. 



0.2 

AVA^ 



(22) 



In the light of the toy model the error in j, for a sample of 
A^ particles at some fiducial radius r, due to Poisson statistics 
scales as aj = a^rl\pN'). So 



;' VorVA^ 
Comparing Eq. ( fT3i and Eq. ( IT9l C, 

(Tj Ao 



= Xq/X implying 



;■ XVN 



(23) 



(24) 



Ao for DM as shown in table above is close to 0.2 which 
explains the scaling relation observed by B2001. 

Distribution of f and its change with redshift 

The distribution of A at any particular redshift is a lojnor- 
mal distribution which can be specified by parameters A and 
<T\. By using Eq. ( I19> we can calculate the distribution of/ 
as shown below 



Pif)df= 



P(X) 



3l/(A,Ao) 



df 



(25) 



z = G.O 
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Fig. 14. — Fraction of counter-rotating matter / as a function of l/<^ = 
v„/(T0 where !'„ is tlie mean tangential velocity and cr^ its dispersion, the 
two parameters of the GSOV toy model described in Section 4. The solid 
line is the relationship / = 1 —Ig(l/Q as predicted by the model. Unlike the / 
vs A plots where there was a large scatter of points about the theoretical curve 
for gas (right column Fig.|5) and a free parameter Aq was required to fit the 
different cases, the data points for the / vs l/i^ plots for both gas and dark 
matter lie on the same predicted curve and they have a very small scatter. 



This predicted formula describes the observed distribution in 
simulations fairly well. The smooth curve in Fig.|5]is the one 
predicted by the Eq. ( I25l l. P{X) does not have a strong redshift 
dependence but / is related to Ao (Eq. ( I19» and Ao increases 
with redshift for gas. This causes the profile of P(f) for gas 
to shift towards smaller values of / at lower redshifts. For 
dark matter Ao is nearly constant so the profile does not show 
any significant change. The fact that < / > for gas decreases 
with decrease of redshift is not only due to increasing ther- 
malization of gas as shown in Section 3.3 but is also due to 
due to two other factors. < / > as given by Eq. ( I16t depends 
on < A >,< A:,- > and < k^ >. < A > increases, while < k,- > 
and < ka > decreases with decrease in redshift. All of them 
act in the same direction to decrease < f >■ 

f as a function ofvo/a^ 

The relation between / and A as given by Eq. ( I19> is an ap- 
proximate one, the actual relation of/ is as given by Eq. Wil . 
Moreover in the plots in Fig.|5lthe gas shows a large scatter. 
So we measure Vg and a^ for each halo and plot / vs Vg/cr- 
These are shown as points in Fig.^J The theoretical predic- 
tion of the toy model as given by Eq. ( I13t is shown as a solid 
line on the same figure. The large scatter which was seen in 
/ vs A plots of gas vanishes and moreover there is no free 
parameter in this relationship. 

Calculating the AMDs by utilizing the toy model 
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Fig. 15. — Comparison of angulai' momentum distributions (P{I) vs / 
plots where / = j/(\/2RyVr)) obtained from simulations with that of profiles 
obtained from the Ga ussi an smearing model (for four different halos same as 
the ones used in Fig. II II . Solid line: profile from simulations Dashed line: 
profile generated by smearing model. The parameters Vo and cr^ adopted for 
generating each of the model profiles is also indicated in each of the panels. 



Fig. 16. — The differential distribution of specific angular momentum for 
various values of a. Both linear and log-log plots are shown. The dashed line 
corresponds to universal profile of Bullock et al. 1 2001 ) with fi= 1.25 and the 
vertical line shows the position of *„„„ for it. a is sensitive to the slope in the 
inner regions, a = 1 marks the transition from a distribution that diverges as 
.V ^ (a < 1) to the one that dips to zero (a> 1). 



The model can also be used to calculate the AMDs. Con- 
sider m(f)dr to be the mass in a cylindrical shell of radius r to 
r+dr. The mean specific angular momentum at this radius is 
given by Vgr. This is smeared by random motion with disper- 
sion (70. So the total mass with radius between r to r+dr and 
specific angular momentum between j to j + dj is given by 
_i — g-(y-i'o»-) /(2(<T0/-) )cijm(r)dr. Integrating this over r we 

get the function m(j) which is the mass of halo with specific 
angular momentum between j to j + dj. 

ft^y 1 

m(j)dj= / -y= e-<>-''°''>'/<2<'"*'''''m(r)c/rc/; (26) 
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For a halo with a given m{r) this integral can be calculated nu- 
merically to give the distribution m{j). We instead follow an 
alternative Monte-Carlo type approach due to its easier imple- 
mentation. For a given halo once the model parameters vv and 
(T are known we calculate the specific angular momentum of 
each of the particles from equation j = (v'o + vv^)^ where vy ^ is 
drawn from a Gaussian distribution of dispersion a. In Fig.llSI 
we have plotted the angular momentum distributions obtained 
by the above procedure along with the AMDs obtained from 
simulations. They are very similar 

5. GENERALIZED PROFILES BASED ON GAMMA DISTRIBUTION 

The various AMDs that have been analyzed in the previ- 
ous sections can be well described by an analytical function 
that depends on just one parameter. The functional form for 
the differential distribution is based on the gamma distribution 
and reads 
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Fig. 17. — The cumulative distribution P(< s) of specific angular momen- 
tum .s = il jtoi for various values of o (left) and /i (right). The shaded area 
separates the region with values of a > 1 form those with a < 1 (unshaded 
area), in both the plots. The a = 1 curve roughly corresponds to curve with 
/J = 1.5. The functional form of [i profiles can extend beyond P(< s) = 1, 
this is shown as dashed lines. The profiles are truncated at s = s,„a\ the point 
where P{< s) = l. For a profiles P(< .v) is always less than 1 and it asymptot- 
ically approaches 1 for large values of s. a and fi measure slightly different 
aspects of the shape, fi curves are a power law with slope of 1 for small s and 
the value of /i is a measure of the point where the bend takes place, a on the 
other hand is sensitive to the slope of the curve for small .v. 



Jd-- 



Jioi 

a 



(28) 



This makes the distribution a one-parameter fit. We choose a 
as the parameter as its effect is easy to understand on physical 
terms. The differential distribution can be integrated to get the 
cumulative distribution as shown below. 



P(j)-- 



1 



JS^a) 



{jr-'e-'ii^ 



(27) 



Since it is normalized it satisfies /^ P{j)dj = 1 . Using the 
fact that /o°° jP(j)dj = J/M = ;„„ = %/2/?,,;,y,.,>A we get 



-Ml- Jo 

where 7 is the Incomplete Gamma function. Writing in terms 
of « = j/jtot and replacing ;,/ we get 



11 



R, = 250.00 kpc c= 10.00 t^= 5.00 kpc ft,= 0.04 _ 
V =124.46 km/s M = 90.05 10'°M„,„ X = 0.0398 



disk 



-M= 



1.05 



Best fit jx and a 
a= 1.41 /i= 2.44 



= 2.0C 




D.Oi 0.10 1,00 10,00 



M= 


2 


^-^^ 


/ a= 1.41 


\ 


/ 


shaded region 
1.3 < a < 1.6 



0.01 0.10 1.00 10.00 



Fig. 18. — The AMD of an exponential disk embedded in an NFW halo. 
Rotation curves are calculated by taking the adiabatic contraction of the halo 
and the flattened geometry of the disk into account. The thick light hne (cyan) 
are the data points corresponding to AMD obtained from the model. A value 
of a = 1.41 seems to fit the model AMD (thin dark; black), /x profiles ai'e 
not very suitable for describing the model profiles. The best fit /.i curve is 
shown as solid line (thin semi dark;red) and also for reference curves with 
^ = 1.05 and /.i = 2.0 are also shown as dashed Unes in the second panel. In 
lower-left panel we plot the cumulative distributions, the solid lines are the 
curves with best fit values of a and ^. The dashed lines here are for curves 
with a = 0.5 and a = 2.0. In lower-right panel the shaded strip corresponds 
to the region where 1.3 < a < 1.6, which is the expected range of values for 
model galaxies with 2 < c < 20 (rj < Tkpc and Ry > lOOkpc). The dots ai'e 
the AMD of the model. 



P{< s) = j(a,as) 



(30) 



Fig. 1161 shows P(s) vs s plots for various values of a in both 
linear and logarithmic plots. For s <C 1/a , P{s) is a power 
law with slope a- I. So for a > 1 P(s) goes to zero as s -^ 0, 
and this gives rise to a dip which is similar to the AMDs of 
observed dwarf galaxies as shown in BBSOl. For a < 1 the 
profiles diverge as i ^ similar to AMDs seen in simulations. 
Fig- El illustrates the effect of a on the shape of P(< s) vs 
s plots, for reference we also illustrate the effect of varying ^ 
(B2001) on the shape of the profiles. 

5.1. What kind of AMD do spiral galaxies have ? 

AMDs of observed dwarf galaxies as illustrated in BBSOl 
are found to have a dip at low s which suggests that they can 
be fit with a profile having a > 1 . The simple case of an expo- 
nential disk rotating in an isothermal halo (p(r) oc 1/r^) with 
flat rotation curve can be calculated analytically and it gives a 
profile with a value of a = 2. To investigate the AMDs of disk 
galaxies in more detail and to find out the range of values of 
a that they satisfy we create a model in which an exponential 
disk of mass fraction fj and scale length r^ is embedded in an 



NFW DM halo of virial rad ius Ry and concentr ation param- 
eter c , similar to analysis bv lMo. Mao. & Whit3 ([Tp98). The 
disk mass fraction is defined as fd = Mdisk/i^hMf/fl,,,). For 
calculating the rotation curve we take into account the adia- 
batic contraction of the halo and also the flattened geometry of 
the disk. So the input parameters for the model are Rv,rd,fd 
and c. Fig.llSlshows a typical AMD obtained by this model. 
The profile is very similar to that observed by BBSOl. It 
can also be seen from the figure that the /i profiles are very 
different from the AMD of model galaxies. For higher values 
of fi (about 2) they can be an approximate fit but still they have 
a distinct core and tail excess. They do not show the charac- 
teristic dip at low /■ as shown by models and also the tail has 
an abrupt truncation whereas the models show a smooth ex- 
tended distribution. So a /i profile even with higher values of 
jj. does not describe the AMD of real galaxies. The general- 
ized profile described in Section 5 provides a good fit to these 
model profiles. For realistic values of the model parameters 
as shown in Fig.^lwe get a ~ 1.4. a is a strong function of 
concentration parameter c but only weakly related to r^ , /?,, 
and fd (keeping other parameters fixed). For lower values of 
fd the fits get more and more accurate, a increases from 1.30 
for c = 2 to 1 .6 for c = 20 (with rj < Ikpc and /?,, > IQQkpc). 
We assume 1.3 < a < 1.6 as the typical range of values ex- 
pected for real galaxies. 

6. TECHNIQUES TO MEASURE THE ANGULAR MOMENTUM 
DISTRIBUTION OF HALOS IN SIMULATIONS 

In this section we first describe two different methods to 
measure the AMDs in halos followed by analysis of AMDs 
obtained by each of these methods. As described earlier the 
velocities given by simulations are of two types. The veloc- 
ity of an SPH particle (gas) refers to the streaming velocity u 
while the velocity of a DM particle samples the actual micro- 
scopic velocity v= u + w, where w is the random motion. One 
method to measure the angular momentum distribution is to 
take the angular momentum of each particle and plot its dis- 
tribution (particle method). This inevitably gives a negative 
tail resulting in a significant fraction of matter with negative 
AM. For gas we find it to be between 4% to 48% with a mean 
of about 15% while for DM it is between 30% to 48% with a 
mean of 41%. This not only makes the interpretation of AMD 
obtained by the particle method difficult but also the compar- 
ison between the AMDs of gas and DM. To make both the 
components come to the same footing we can either broaden 
the velocities of gas particles or smooth the velocities of both 
gas and DM particles over a fixed number of neighbors. How- 
ever, even smoothening the velocities over 1000 neighbors is 
not enough to suppress the negative angular momentum tail. 

B2001 obtain the AMD by dividing the halo into cells and 
then calculating the angular momentum by averaging over all 
particles in the cell (cell method). These cells are then used 
to plot the cumulative distribution. In the particle method de- 
scribed above it is possible to derive both the differential and 
cumulative profiles while in the cell method only latter is pos- 
sible owing to the small number of cells, typically around 60. 
Furthermore, the cells are assumed to have a spherically sym- 
metric geometry. Each cell covers a full lir range in cj) and 
they span the range of (r/Ry,sin9) from (0,0) to (1,1). The 
radial shells are spaced such that each of them contains the 
same number of particles. The shells are then divided into 
3 azimuthal cells of equal volume between sin9 = and 1 . 
Positions with same rsin9 above and below the plane belong 
to same cell. For our calculations we divide the halo into 
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Fig. 19. — The cumulative angular momentum distribution for 6 different 
halos obtained by the particle method (dark solid line). For each halo, distri- 
butions for both gas and dark matter are shown. The best fit a and fj, values 
are also labeled . Light solid line (orange): best fit fj, profile constant A fit. 
Dashed line (magenta): best fit a profile. Up down aiTows mark the point 
with maximum value of .v. The profiles have a smooth truncation the slope in 
outer parts gradually goes to zero. i,„av is higher than that obtained by cell 
method. 



Fig. 20. — The cumulative angular momentum distribution for 6 different 
halos obtained by symmetrical cell method (dark solid line). For each halo, 
distributions for both gas and dark matter are shown. The best fit a and fj, val- 
ues are also labeled. Light solid line (orange): best fit fj, profile constant A fit. 
Dotted line (red): best fit fi profile constant jmax fit. Dashed line (magenta): 
best fit a profile. Up down arrows mark the point with maximum value of i. 
The constant jmo, fits always pass through the point with maximum s while 
the constant A and a profiles in general do not. The profiles tend to have an 
abrupt truncation, which means that at i = s,„ai the slope is significant. 



60 cells with approximately constant number of particles in 
each of them. We call this method symmetrical cell method. 

The symmetrical cell method is more effective in reducing the 
negative AM material compared to particle method. However 
the resulting AMD might be biased if the system is not axis 
symmetric. We therefore repeat the analysis with a modified 
scheme which is free from any inherent symmetry. We divide 
the halo into 4 radial shells with the nth shell containing 2n^ 
particle (2,8, 18 and 32). Each shell is divided into n polar 
zones denoted by I. Each / region is divided into 2/+ 1 az- 
imuthal zones. Finally regions above and below the z plane 
belong to different cells. The radial and polar divisions are 
done such that each cell contains a constant number of parti- 
cles. We refer to this method as normal cell method. 

6. 1 . Analysis by cell method and Bullock profiles 

In B2001 cumulative AMDs were calculated by the sym- 
metrical cell method for DM halos and a universal profile with 
shape parameter ji was shown to fit the data. The distribution 
of ^ was found to satisfy a log-normal distribution with < 
logiaip,- 1) >= -0.6 and a = 0.4. Our simulation s reproduce 
these findings as shown in Fig.|23](second row). Chen & Jind 
(|^002, henceforth CJ02) have also obtained similar results ex- 
cept the fraction of cells with negative AM is higher than in 
B2001. B2001 found that 5 % of halos have f > 0.1 while 
CJ02 found it to be 40%. IChen. Jing. & Yoshikaw (2003 
henceforth CJ03) have also analyzed the AMDs for both gas 
and DM components. They found that /x is higher for gas than 
that for DM and also that / is lower for gas. Results of our 
simulations agree quite well with the findings of CJ02 and 
CJ03. We find that for DM / > 0. 1 for 33% of halo while for 
gas / > 0. 1 for 9% of halos. < / > for DM is -- 9% while for 
gas it is ^3%. We also find ji to be higher for gas compared 




Fig. 2 1 . — The cumulative angular momentum distribution for 6 different 
halos obtained by normal cell method free from any symmetry restrictions 
(dark solid line). For each halo, distributions for both gas and dark matter are 
shown. The best fit a and /.t values are also labeled. Light solid line (orange): 
best fit /t profile constant A fit. Dotted line (red): best fit ^ profile constant 
jmax fit. Dashed line (magenta): best fit a profile. Up down arrows mark the 
point with maximum value of s. The profiles end at higher values of .v and 
they also tend to have smoother truncation compared to profiles obtained by 
the symmetrical cell method. They are more similar to the profiles obtained 
by the particle method. 



13 



.000 



0.100 



0.010 



IJ.^ 1.52 a=1.09 
Particle/GAS 




/i,= 1.35 a = 0.94 

Particle/GAS 
halo-2 



V 0.100 




,„=1.35 fJ.t=^.^jl a = 0.99 

ill/GAS 
lolo- 1 




jii= 1.09 c< = 0.75 
■-'■■— -JParticle/DM 
;nolo-l 



j„„=1.74 Mx=1.52 0=1.02, 

;ell/GAS 
ialo-2 




/i= 1.12 a = 0.76 
■— — HParticle/DM 
;nalo-2 




„„=1.36 |Uj=1.33 a = 0.90 
ell/DM 

;halo-2 




.1 1 a = 0.63 



/i.= 1.09 a = 0.72 

Particle/GAS 
holo-3 



0.100 



0.010 




- -1E 



,„=1.27 Aij=1.16 a = 0.69 

iell/GAS 
ialo-3 




/i= 1.08 a = 0.74 
-|=article/DM 
;nalo-3 




/fjmox=l-17 Mi=1.18 a = 0.79 

iEell/DM 

;holo — 3 




^l= 1.11 a = 0.75 
Particle/GAS 



V 0.100 




IJ.= 1.04 a = 0.63 
Particle/GAS 
tialo-5 



0.010 



j„„=1.51 /!j,j=1.40 a = a.i 

iell/GAS 
ialo-4 



/i= 1.05 a = 0.66 
-pQrticle/DM 



,= 1.22 |Ui=1-18 a = 0.63 
:Eell/DM 
:tialo-4 





imo,= 1^1+ Mi='-08 a = 0.60 

iell/GAS 

ialo-5 




/i= 1.04 a = 0.63 
-HPorticle/DM 
;nalo-5 




„„=1.13 /j,j=1.15 a = 0.65 
ell/DM 
itiolo — 5 




/j,= 1.09 = 0.72 

Particle/GAS 
tialo-6 



V 0.100 



0.010 




ftj„„=1.16 /.tj=1.17 a = 0.i 
^Eell/DM 
;halo — 6 




0.001 

0.001 0.010 0.100 1.000 10.000,001 0.010 0.100 1,000 10.000.001 0.010 0.100 1,000 10.000.001 0.010 0.100 1.000 10,000 



Fig. 22. — The cumulative angular momentum distribution for 6 different halos (dark solid line). Each row corresponds to a particular halo. Profiles of each 
halo both for gas and dark matter by particle and as well as symmetrical cell method are shown. Also are labeled the best fit a and /t values. Light solid line 
(orange): best fit [i, profile constant A fit. Dotted line (red): best fit /i profile constant jmax fit- Dashed line (magenta): best fit a profile. The subscript to [i, denotes 
if it is a constant A or a constant j,nax fit. For particle method the ^ profiles are generated by constant A fits only. 
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Fig. 23. — Disti'ibution of shape parameter fi for angular momentum 
distributions obtained by various different methods and fitting techniques. 
Median, mean and standard deviation of logip, - 1 ) is also indicated for each 
case, and below it coiresponding /i value is shown. For gas the constant A 
fits have lower mean and higher width compared to constant jma_\- fits. On the 
other hand for DM the mean increases. Same trend is observed in analysis 
with normal cells. In comparison with symmetrical cells the normal cells 
give a lower mean and width for the fj. distributions. The particle method 
gives even lower mean but the widths are maximum. 



to DM (second panel Fig.l23> by using the same methods for 
analysis as done by above authors. 

However, we observed that the /i values obtained from fit- 
ting M(< j) vs j data are sensitive to the fitting procedure 
and error bars used. In fact two techniques can be used to fit 
fi profiles each giving slightly different results. The analytic 
function for the /i profiles is given by 

M(< /■) ni 

^Tr^=Pi<j) = —r- (31) 

My JO + J 

jo can be written in terms of A and /i as 

V2VyRA 



J0-- 



Jtot 



(32) 



b(p) b(p) 

where bin) = -jjl ln{\ -/i"')- 1. This makes the jjl profiles a 
one parameter fit for a halo with a given A. 



Pi<j)-- 



P«s) = 



N 



;o(/i,A) + ; 



(33) 



(34) 



l/b(^i) + s 

We call this the constant A fits. Alternatively the ^ profiles 
have an implicit maximum specific angular momentum j,„„^ = 
jo/ifJ.- !)■ Writing in terms of j /],„„_, 

A* J / Jmax 



P«j)-- 



(fJ.-l) + j/j„ 



(35) 



For a given halo in simulations /„„„ is the specific AM of 
the cell having maximum specific AM. This is again a one 
parameter fit which we call constant j^ax fits. 

In B2001 it was shown that for DM the error bars on values 
of j for DM can be approximated by cry = Q.2/{X^/(N)). The 



Fig. 24. — Distribution of shape parameter a for angular momentum dis- 
tributions obtained by the particle and by the cell method. Median, mean and 
standard deviation of log(a) is also indicated for each case, and below it cor- 
responding a value is shown. The vertical line corresponds to a = 1.3. The 
normal cell method gives lower values of mean and width for both gas and 
DM components compared to symmetrical cell. The same trend is observed 
for particle method except that the mean value for gas increases. 
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Fig. 25. — ^oas - 1 vs iMoM — 1 for various different techniques. Median 
and mean of (fiQas ~ 1)/(Mdm ~ 1) ^re also labeled on each of the plots. The 
straight line shown in the plot corresponds to the relation fiQas = f^DM ■ Parti- 
cle method and constant j„,ax fits have ficas > fJ-DM but constant A fits do not 
show any significant bias. 



velocity dispersions are in general not homogeneous through- 
out the halo. CJ03 therefore base their analysis on error bars 
estimated by aj = Vc(r)r/\/(N). However, gas particles have a 
much smaller velocity dispersion. We therefore use a scheme 
in which the error are estimated by aj = jstd-dev/VN, where 

j^d-dev ~< 0~ < i >)^ > is the standard deviation calculated 
over the particles in the cell. For the cumulative distribution 
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Fig. 26. — OLGm vs ciDM for various different techniques. Median ,mean 
and standard deviation of (a-Q^sl'^DM) are also labeled on each of the plots. 
The straight line shown in the plot corresponds to the relation O-cas = cidm ■ 
In particle method we see uq^^ > aoM but cell method does not show such a 
bias. The two vertical and horizontal Hues correspond to a = 1 .3 and there are 
very few halos that lie above these lines (except for the middle panel which 
has a few). 



Fig. 27. — Comparison of parameters fans and foM for various different 
techniques. < fgas > and < foM > are also shown. The distribution for 
normal cell method is very similar to that of particle method. In symmetrical 
cell method / is in general much lower and for gas a significant number of 
the halos have / = 0. 




Fig. 28. — Effect of broadening the velocities on the distribution of angular 
momentum for two different halos. After broadening the angular momentum 
distribution of gas is very similar to that of dark matter. 



with CJ03 (Fig. |23]and Fig. 1251. second row). The distribu- 
tion of logifi- 1) can be roughly fit by a Gaussian and the fit 
parameters for DM are similar to those in B2001. For con- 
stant A fits however, we find that figas is nearly same as udm- 
The reason for the discrepancy is that the constant j,„ax fits are 
constrained by construction to satisfy P{< j) = 1 at y = j,„ax, 
so the error in the value of j„i„x is not taken into account in 
the fitting procedure. Moreover the fits are not constrained to 
satisfy Eq. ( I32> . For constant A fits and also the a fits this is 
not the case, the curves do not necessarily truncate at j = j,„ax'- 
they may extend to j > j„,ax or may already stop at j < j^ax 
(Fig.l20>. There is a systematic trend such that for cases where 
/i is high (in constant j,„ax fits), the constant A fits truncate at 
j > jmax, resulting in a lower value of /i. The DM does not 
have high values of /i so it is relatively unaffected while gas 
has relatively high values of /i and is significantly affected. 
In other words, the effect that Ugas > fJ-oM is diminished to a 
large extent in the constant A fits. The distributions of agas 
and aoM also show only a mild bias, (Fig. 1241 second row ). 
About 20% of halos had a > 1 .3 but out of these for 10% of 
halos the fits were poor so they were rejected. After correc- 
tion only 10% of halos have a > 1.3 for gas, while for DM 
the percentage is about 5% (Fig.l26>. 

Changing the geometry of the cells changes the distribu- 
tions of /i and a slightly. The peak position and width of the 
distributions are both reduced. A smaller percentage of halos 
have a > 1.3. We address the reason for these results in the 
next section. 



^(< j) vs j, <Jj gives the error along j-axis. We estimate the er- 
ror along f< axis by ap^ = {dP^/dj)aj. To calculate dP^/dj 
we first fit the data with a /i or a profile and then use this 
value of fioi a to calculate dP^/d j. AMDs obtained by both 
particle and cell methods along with corresponding fi and a 
fits are shown in Fig.[^ Fig.|20| Fig. 12 Hand Fig. 1221 
We find that for constant j^ax fits iigas > I^-dm in agreement 



6.2. Analysis by particle method 

Is broadening technique suitable for comparing AMD of gas 
and DM? 

It was shown in vB2002 that the AM profiles of gas and DM 
are remarkably similar, if not identical, when the velocities 
of gas were broadened by their microscopic thermal motion 
('Fig. l28t . However, if the shape of the profiles after broaden- 
ing is determined primarily by the value of a and is thus insen- 
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Fig. 29. — Angular momentum distribution of gas of six different lia- 
los before broadening (top panel) and after broadening their velocities with 
a = 160.1 kms"' (isothermal broadening) (lower panel). After broadening 
all the different distributions are similar so broadening is not a very reliable 
technique for making comparison between the distributions of gas and DM. 



sitive to the original shape, then different profiles can be made 
to look similar by broadening with same a. We demonstrate 
this in Fig. 1291 where we broaden the streaming velocities of 
6 halos, each with a unique AMD, with same velocity disper- 
sion cr= 160.1kms"'. Although the initial profiles were quite 
different, the profiles after broadening are very similar We 
conclude that thermal broadening masks out the uniqueness 
of un-broadened AMDs and hence is not a suitable technique 
for making comparisons. 

Angular momentum distributions after smoothening 

Rather than broadening we consider smoothening of the 
DM velocities as the superior procedure, i.e. we extract 
streaming velocity u from the given velocity v and then com- 
pare their AMD with that of gas. Increasing the numbers of 
neighbors N„gb, over which the smoothing is applied, will de- 
crease / but then the region over which averaging is done 
eventually becomes too large and any significant informa- 
tion is lost. Furthermore, particles close to the rotation axis 
suffer an apparent loss of AM due to a geometrical effect: 
< v^ > and < Vy > in the inner region are zero for a symmet- 
ric rotating system aligned along z axis. Alternatively we may 
smooth the AM instead ( the halo should be centered before 
smoothening). In this case we encounter the following prob- 
lem: Particles close to the axis and in a conical region around 
it get enhanced in AM while those along the equator suffer 
a loss in AM due to the existence of a strong density gradi- 
ent. A spherical volume around such a typical particle has 
more neighbors towards the symmetry axis with lower AM 
than away from it, so the smoothed AM which has a radial de- 
pendence given by j{r) = Vgr, is lower close to the axis. This 
effect was not prominent when the smoothening was done on 
velocities because the velocity does not exhibit a strong radial 
dependence (v(r) ^ v„). This suggests that a better method 
would be to smooth the velocities in cylindrical co-ordinates. 
For our analysis here we use the technique of smoothening the 
angular momentum only and we choose the number of neigh- 
bors to be 400. Instead of an SPH type kernel we use a simple 
step function, which is equivalent to taking the mean over the 
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Fig. 30. — Effect of smoothening o n the angular momentum distributions 
of two different halos (same as in Fig. 1281 . The angular momentum of both 
gas and dark matter is smoothed by taking mean over 400 neighbors. For 
some halos after smoothening the gas and DM profiles are similar while for 
others they are different, one of the reasons for this is that the A and / are not 
same for gas and DM. 




Fig. 31. — Effect o f sm oothening on the P(s) vs s plots of two different 
halos (same as in Fig. I28i . The angular momentum of both gas and dark 
matter is smoothed by taking mean over 400 neighbors. For the halo in the 
lower panel, the gas and DM profiles after smoothening are similar but for 
the other halo they are different. 



neighbors. This makes it easier to make comparisons with the 
cell method and is also more effective in reducing /. 

In(Fig. I30> the effect of smoothening on the distributions 
of angular momentum is shown. For some halos the profiles 
of gas and DM are similar after smoothening and for some 
they are different. This kind of comparison is not very useful 
because the shape of P(l) vs I plots depends on A and /. So 
a good agreement may merely reflect that the A and / are 
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Fig. 32. — A comparison of fi and a fits appHed to angular momentum 
distributions of the gas obtained by the particle method. The vertical line in 
plots on the right indicate the truncation point s„,a.i of fi profiles. The solid 
lines are the data obtained from simulation and the dashed lines are the best 
fit a (left) and fj. (right) profiles. In outer parts the actual profiles are steeper 
than 1/.V" (the expected asymptotic form of /t profiles). 



similar for gas and DM. Moreover A and / are in general not 
same for gas and DM, this makes the interpretations of these 
plots even more difficult. 

So to make a comparison first we need to take out the de- 
pendence on A and / which can be done by plotting P(s) vs 
s for the positive tails of gas and DM and then comparing the 
best fit values of /i or a. This is shown in Fig. 13 II 

Analysis of smoothened profiles 

We bin the smoothed profiles such that each bin contains 
Nhaio/^OO particles. For calculating the error bars we use 
the same technique as described earlier in the case of the 
cell method, with the exception that we assume aj/j =< ^ > 

/VN. For DM < ( >= 5 and with A^ = 400 this gives a value 
of 0.25. For gas < (^ >= 1 and this gives a value of 0.04. 
The smoothed profiles can be fit by both Bullock /i profiles 
(only constant A fits are used) and the generalized a profiles 
Fig- 13 P(i) vs s differential plots in Fig. 1321 show that that 
the Bullock profiles are shallower for large s, i.e. the actual 
AMDs are steeper than l/s^ in outer parts. The particle pro- 
files do not truncate abruptly like in the cell method, the slope 
gradually goes to zero. The a profiles provide better fit to the 
particle profiles in outer regions. 

The /i and a values obtained by particle method are also 
larger for gas as compared to that of DM (Fig.l25land Fig.l26t. 
logip,- 1) and log(a) distributions are roughly Gaussian. The 
mean, median and standard deviation are shown in Fig. |23] 
and Fig. 123 a > 1.3 for about 2% of gas halos while none 
of the DM halos have a > 1.3. Even after smoothening with 
same number of neighbors in particle method / on average is 
greater for DM (13%) than for gas (8%). Similar effect is also 
seen in cell method where the cell size used is same for both 
DM and gas Fig. 1771 

Effect of spherical and symmetric cell geometry 



A comparison of results of particle method with that of 
symmetrical cell method show that /i (constant A fits) and 
a are higher for symmetrical cell method, ii specifically is 
higher because s„,ax to which it is directly related is about 
1 .5 - 2.0 times higher in particle method than in symmetrical 
cell method (Fig.ll9IFig.l20ll although the number of neigh- 
bors in particle method is approximately the same as the num- 
ber of particles in each cell. / is also significantly lower in 
symmetrical cell method than in particle method (Fig. I27t . 
Since smoothening is a local effect whereas cell averaging is 
not due to its spherical and symmetric construction. So a blob 
of negative AM or higher AM material will retain its charac- 
ter when smoothed but will be averaged out in cell method 
over regions that are not necessarily local. Both methods will 
deliver identical results only in the case where the system is 
symmetrical. Any deviation from symmetry, in the symmetri- 
cal cell method has the effect of lowering the low AM material 
by mixing it with high AM one and vice versa. So both / and 
Smax are lower in symmetrical cell method, a does not depend 
upon s,„ax but is only affected by the lowering of low AM ma- 
terial and this effect is not very strong so it shows only a slight 
change. 

These arguments also imply that the results of normal cell 
method should be very similar to particle method. This is in- 
deed the case as can be seen in Fig.l27lwhere distribution of/ 
for normal cell method is very similar to particle method. The 
distribution of a and /i parameter for particle method is also 
closer to normal cell method than the symmetrical cell method 
as shown in Fig.|26l In spite of the similarities there are subtle 
differences specially for gas whose < log(a) > is much lower 
for the case of normal cell method. This discrepancy may be 
because the number of particles in the normal cell is not fixed 
for all halos like in particle method for which the number is 
400. Compared to normal cell method the particle method has 
slightly lower values of /i because their s„ax is slightly higher 
and this is because AM is a monotonically increasing function 
of radius and for a cell of finite radial thickness the average 
AM of the cell will always be less than the maximum AM of 
particles in it. 

Measuring the spatial asymmetry of angular momentum 
distribution within a halo 

The above discussion suggests that both gas and DM should 
have significant asymmetry. We verify this as follows: We 
measure the asymmetry by dividing the halo into cells as de- 
scribed earlier but this time each radial shell is also divided 
into 6 azimuthal zones and 3 9 zones. Then we measure the 
z component of angular momentum vector jir) and jir') for 
a pair of cells with opposite parity situated at r and r' where 
r' = -r. Symmetry of j distribution within a halo is given by 

„ Jj(r) + j(r')\^-\j{r)-j{r')\^ ^ 

O ;■ = < 7— T—TT^ TT— 7— TTT >all cells 



' ^ \i{r) + j{r')\^+\j{r)-j{r')\^ 
2j(r)j(r') 



= < 



jHr) + jHr') 



>all cells 



(36) 



Sj can vary from -1 (perfectly anti-symmetric system) to 1 
(perfectly symmetric system). An asymmetry as high as 
j{r) = 2j(r') corresponds to Sj of only 0.8. < Sj > ,the mean 
over all halos, is found to be 0.78 for DM and 0.83 for gas. 
So in halos both gas and DM have significant asymmetry as 
expected. DM in fact is more asymmetrically distributed than 
gas. 
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Table 4. Comparison of x for fi and a fits 



Cell method " Particle method 



DM GAS DM 



GAS 



median(xi/xp 0.78 1.31 
0.25 0.27 



median( —2 — tL ) 



1.47 0.92 
0.39 -0.08 



"Comparison shown for symmetrical cell method only 



Table 5. Pearson correlation coefficient of a with other halo parameters 



M, 



Sj 



syiii 



par 



-0.12 

-0.16 



0.27 
0.53 



-0.14 
-0.25 



-0.04 
0.27 



0.51 
0.76 



0.86 
0.80 



0.69 
1.00 



Note. — The superscript on a denotes the method employed to calcu- 
late the AMD e.g particle method , symmetrical cell method or normal cell 
method (free from symmetry restriction). 



6.3. Comparison of quality of a and fifits 

It might sound strange that two different functional forms 
are being used to describe the same data. But it is not very 
surprising considering the fact that the profiles cannot be per- 
fectly described by a functional form, they have small devi- 
ations and this gives enough room to both p or a profiles 
(which after all are not very different) to be used to fit them 
with equally good fits. The profiles from the cell method 
(symmetrical) tend to have an abrupt truncation at j = j,„ax 
for some cases ( probably due to the symmetry effects dis- 
cussed in previous section and also because the number of 
data points used to sample the profile, which is the total num- 
ber of cells is small, typically around 60 ). So /i profiles which 
by design have an abrupt truncation at j = j,„ax perform better 
in the outer regions for these cases, but for the same reason 
they are not able to fit the particle profiles in the outer regions 
which have a smooth extended tail. The a profiles fare better 
for these. So if we neglect the outer parts P(< s) > 0.95 then 
both the profiles provide a satisfactory fit to the distributions 
from both methods. To get an estimate of the quality of the 
fits, we measure x^ for data points with P{< s) < 0.95 for 
both profiles, x^ is defined as x^ = Eti((3'f '""3'"'"'"')/^;.)^ 
for each of the fits. To make a comparison we calculate the 
median of (Xa -X^)/(XaXp) and median of Xa/x^ over all 
the halos. If the x^ for both of them are equal then the for- 
mer quantity is close to zero. A negative value of -1 implies 
X?i = 2.5Xq- If the median of Xq/Xu is close to 1 then both fits 
are good for equal percentage of halos. If more than 50% of 
halos have Xa < Xp then the above quantity will be less than 
1 and vice versa. In Table Oand Table |5]the values of these 
quantities are listed for various methods. Form the table we 
can see that both the profiles are equally good at describing 
the AMDs obtained by symmetrical cell and particle methods 
(except for the outer 5% ). 

7. CONCLUSION & DISCUSSIONS 

We have presented here results from the non-radiative hy- 
drodynamical simulation of 41 high resolution halos whose 
masses were selected to span the range from dwarf to bright 
galaxies. Our investigation mainly focused on the angular mo- 



mentum properties of the halos and whether there are system- 
atic differences between gas and dark matter. Our findings 
can be summarized as follows: 

1 . We investigated some of the global angular momentum 
properties like spin parameter A, fraction of negative 
angular momentum / and misalignment angle 9. We 
find the spin parameter of gas to be on average larger 
than that of DM and this effect is systematically more 
pronounced at lower redshifts. At z = \gas/^DM ^ 1 -4, 
which is in agreement with the result reported by CJ02, 
but is not in agreement with vB2002 who do not find 
any such bias. The mean of the misalignment angle 9 is 
20° which is again in agreement with CJ02 who get a 
value 23.5° but is less than the value of 36.2° obtained 
by vB2002. Both these discrepancies could be due to 
the inclusion of a large number of low resolution halos 
in the analysis of vB2002. The counter-rotating frac- 
tion / is anti-correlated with A and for gas / decreases 
with decrease of redshift, an effect that can be explained 
by the increasing level of thermalization at lower red- 
shifts. Other than this there is little evolution of other 
properties with redshift. 

2. We find that the fraction of material with negative an- 
gular momentum can be described by the equation / = 
1 -Ig{\/\o), with I{x) being a Gaussian integral. To un- 
derstand these effects we developed a toy model, where 
we introduce an ordered velocity Vo which is smeared 
by means of Gaussian random motion with dispersion 
a. This model reproduces the / = /(A) correlation and 
suggests an actual relation of / = 1-/^(1/0 where 
C = cr/vo- This relation is demonstrated to be in excel- 
lent agreement with the results from simulations. We 
also see that the gas gets more and more thermalized at 
lower redshift resulting in cr/V,, to be smaller at lower 
redshifts, whereas A and concentration c increases at 
lower redshifts. All these effects are contributing to 
the decrease of the amount of gas with negative angular 
momentum. The model also reproduces well the shape 
of AMD observed in simulations. 

3. We study the distribution of angular momentum in de- 
tail, and compare and contrast various different tech- 
niques used to derive angular momentum distributions. 
We first use the cell method as proposed by B2001 and 
reproduce the result of CJ03 that jjl for gas is greater 
than that of DM. We find that the results are sensitive 
to the very details of the method employed for fitting. 
In particular the effect of gas having higher ji is dimin- 
ished to a large extent if fits are performed assuming 
a constant A rather than a constant y'max- According to 
CJ03 fjL for gas is comfortably in the range required by 
observation of disk galaxies which is (p. > 1.75). We 
also find that about 30% of halos have /i > 1.75. 

4. By comparing the AMD found in the simulations with 
those of exponential disks we conclude that merely hav- 
ing /i > 1.75 is not a sufficient condition to match the 
angular momentum profiles of observed disk galaxies. 
We find that a generalized profile, based on gamma dis- 
tribution, with a single parameter a can be used to fit 
the AMD of model galaxies (exponential disks embed- 
ded in NFW halos) as well as AMD of gas and DM in 
simulations, a > 1.3 seems to resemble the profiles of 
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dwarf galaxies shown in BBSOl, a condition that is only 
obeyed by a small minority of halos (less than 10%). 
For fits based on the (symmetric) cell method only 10% 
of the halos have gas with a > 1.3 (after rejecting ha- 
los with bad fits). For particle method and normal cell 
method the percentage is even lower, about 2% of ha- 
los have a > 1.3. We find the particle method to be 
more robust, with data that is less noisy and also free 
from any artifical non local averaging or any symme- 
try assumptions. The profile in particle method do not 
have abrupt truncation like in cell method (symmetri- 
cal) and are more extended. This may have important 
implications for the truncation radius and extent of gas 
in real disks: In semi analytical models the distribu- 
tion and extent of cold gas depends upon the AMD 
used for the models. In van den Bosch ( 2001 )(Fig-9 
there) the models have an AMD with fj, between 1.6 
and 1.9, and the disks are predicted to have a sharp 
truncation {Rgas/RHi ^ 1) of cold gas which is in dis- 
agreement with observed distribution of HI in galaxies 
{Rgas/RHi > 1-5). Here Rgas is the maximum extent of 
gas with non zero surface density and Rhi is the radius 
with surface density of 1 Mq/jc"^. 



5. For the particle method, there is a significant fraction of 
counter-rotating matter which has been excluded in cal- 
culating a. The final AMD (and therefore the predicted 
structure of the model galaxy) will depend strongly 
upon how this material eventually gets mixed up with 
the remaining portion of the halo, during the assembly 
of the galaxy. We saw that the counter-rotating mat- 
ter undergoes extensive mixing during the history of its 
formation in the hierarchical framework. But the sit- 
uation we have here is slightly different as it concerns 
the fate of the counter-rotating gas in an isolated halo 
when subjected to collapse. In the absence of any infor- 
mation it is best not to assume any preferential mixing 
with low, or high AM material or also equal AM mate- 
rial in which case it might form a non rotating bulge as 
prescribed in vB2002. The most plausible prescription 
seems to be that of random mixing: During collapse the 
counter-rotating matter will follow a slightly different 



trajectory than normal matter around it and in the pro- 
cess it will get mixed by shocks with matter of various 
different AM along its path. This process does not seem 
to have any preference. If we assume that the process 
of mixing will not change the distribution of a signif- 
icantly, we can conclude that less than 10% of halos 
have a > 1.3. In other words only the absolute minor- 
ity of halos have AM distributions that resembles that 
of an observed (dwarf or LSB) disk galaxy. In the ab- 
sence of any significant correlation of a with mass M,, 
or for that matter with any other halo parameter, this 
leads us to the general conclusion that a typical halo in 
ACDM simulation has far too much low angular mo- 
mentum material to account for majority galaxies fea- 
turing a dominated disk component. 

6. If the angular momentum of each element is conserved 
then low angular momentum material has to be pref- 
erentially discarded during the process of galaxy for- 
mation, this might be effected by means of supernova 
feedback during star formation that drives out gas from 
small halos that come randomly from all direction and 
contribute mainly to the low angular momentum ma- 
terial as shown by Mailer & Dekel (2002). The AMD 
of bright galaxies in Mailer & Dekel (2002), which are 
effected little by feedback , seems to have a ^ 1 ( no 
dip in AMD for small s), consistent with AM profile of 
a galaxy with a bulge which accounts for the low AM 
material. It is no way similar to exponential disk with 
flat rotation curves as claimed there, because that would 
imply an a = 2 and will be reflected as a prominent dip. 
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